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Abstract 

In this paper we address the topic of inelastic electron scattering in meso- 
scopic quantum transport. For systems where only elastic scattering is 
present, Landauer theory provides an adequate description of transport that 
relates the electronic current to single-particle transmission and reflection 
probabilities. A formalism proposed recently by Bonca and Trugman facil- 
itates the calculation of the one-electron transmission and reflection prob- 
abilities for inelastic processes in mesoscopic conductors connected to one- 
dimensional ideal leads. Building on their work, we have developed a self- 
consistent procedure for the evaluation of the non-equilibrium electron distri- 
butions in ideal leads connecting such mesoscopic conductors to electron reser- 
voirs at finite temperatures and voltages. We evaluate the net electronic cur- 
rent flowing through the mesoscopic device by utilizing these non-equilibrium 
distributions. Our approach is a generalization of Landauer theory that takes 
account of the Pauli exclusion principle for the various competing elastic and 
inelastic processes while satisfying the requirement of particle conservation. 
As an application we examine the influence of elastic and inelastic scattering 
on conduction through a two site molecular wire with longitudinal phonons 
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using the Su-Schrieffer-Heeger model of electron-phonon coupling. 
PACS: 73.50.-li, 73.23.-b, 73.61.Pli 
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I. INTRODUCTION 



The role of inelastic scattering in electron transport through mesoscopic systems is a 
topic of current theoretical and experimental interest. Recent experiments have revealed 
the importance of the electron-phonon interaction in transport measurements performed on 
certain mesoscopic systems. In particular inelastic scattering effects have been observed 
directly in STM measurements of the differential conductance of molecules adsorbed on 
metallic substrates.B^ Theoretically the effects of inelastic scattering on electron transport 
through mesoscopic semiconductor devices have been investigated by a variety of methods 
ranging from Green's function techniqueJl^ to the Fermi golden rule.B Some theoretical 
models have also been proposed to elucidate the effects of molecular vibrations on electron 
tunneling through molecules contacted by an STM.Hl^ Recent work has provided a method 
for calculating the electronic current at finite temperatures for electrons tunneling through 
a one-dimensional conductor in the presence of phonons.0 These approaches have yielded 
valuable insights into the behavior of specific systems. 

In the absence of inelastic scattering and electron correlations, Landauer theoryS pro- 
vides a general framework for calculations of the electronic current through mesoscopic 
conductors that are coupled to ideal single- or multi-channel quasi-one-dimensional leads. It 
relates the electronic current to the transmission probability for an electron incident from the 
source lead to scatter elastically through the conductor and into the drain. The transmission 
probability is found by solving the single-electron quantum scattering problem. 

When inelastic processes such as phonon emission and absorption are considered, electron 
scattering becomes a many-body phenomenon involving electrons and the various excited 
phonon states of the system. A non-perturbative approach has been proposed by Bonca 
and Trugman for treating this scattering problem for a system consisting of a mesoscopic 
conductor (which supports the phonon modes) coupled to two single-channel ideal leads 
that act as the electron source and drain.0 Their method approximates the many-body 
problem by a multi-channel single-electron scattering problem which can be solved exactly. 



3 



Here each channel corresponds to a different vibrational quantum state of the mesoscopic 
conductor. Using their approach it is possible to determine the transmission and reflection 
probabilities for all inelastic and elastic scattering events that an electron may suffer in the 
conductor as it goes from the source to the drain. This method has been applied to study 
the effects of phonons on electron transmission through one- dimensional conductors where 
the electron-phonon coupling has been modelled using the Holstein HamiltonianBffl and 
the Su-Schrieffer-Heeger (SSH) Hamiltonianli^lizl. 

The present article is complementary to the above theoretical work. Its purpose is 
to generalize the Landauer theory of electrical conduction to mesoscopic systems with 
electron-phonon scattering, assuming that the latter is described accurately within the 
multi-channel single-electron scattering approach of Bonca and Trugmanlli. I.e., given the 
Bonca-TrugmanS solution of the inelastic scattering problem, our objective is to calculate 
the electric current by appropriately generalizing Landauer theory. 

When transport in a many-electron system is described in terms of one-electron scattering 
processes, the Pauli exclusion principle needs to be considered since in general electrons 
involved in different transitions may compete to occupy the same final state. In treatments 
of transport based on the Boltzmann equation, if the electron scattering is elastic, this 
competition has no effect because of a mutual cancellation of terms.lii However such a 
cancellation does not occur when scattering is inelastic.0 In Landauer theory scattering is 
assumed to be elastic and the electronic current is carried by electrons which occupy single- 
particle scattering states that extend across the mesoscopic conductor from the source to the 
drain. These scattering states are orthogonal to each other.0 Because of this orthogonality, if 
an electron populates one of these scattering states it does not compete with any of the other 
filled scattering states. Thus in Landauer theory (as in the elastic case of the Boltzmann 
equation) the Pauli exclusion principle for the final states in the scattering processes does 
not play a role in the determination of the electronic current. 

In our generalization of Landauer theory that includes inelastic scattering, this is no 
longer true: Different scattering processes can send electrons to the same outgoing state and 



the Pauli principle plays an important role in determining the occupations of the outgoing 
channels. In the Boltzman treatment of inelastic scattering the Pauli principle is applied 
to the non-equilibrium one-electron distribution function for the system since transport is 
an inherently non-equilibrium phenomenon. Here we apply an analogous principle: In our 
extension of Landauer theory the Pauli exclusion principle is applied to the non-equilibrium 
electron distributions in the outgoing channels of the scattering process and these distribu- 
tions must be determined self-consistently. Since some processes may be excluded because 
of the Pauli principle, overall particle conservation is another consideration that must be 
taken into account when trying to generalize Landauer theory to include inelastic scattering. 

We have developed a practical method for determining these non-equilibrium electron 
distributions for a class of mesoscopic systems. These systems are similar to those considered 
by Bonca and TrugmanEl in that they consists of a mesoscopic conductor which is attached 
to two single mode ideal leads which act as a source and drain, and phonons are assumed to 
be present only on the conductor and not on the ideal leads. We consider scattering processes 
which change the conductor's phonon state or leave it unchanged. Each of these processes 
contributes to the non-equilibrium electron distribution. We determine each contribution 
self-consistently (i.e., in the presence of the other processes that can occur) in terms of 
the transmission and reflection probabilities for the given process and the Fermi functions 
that describe the incoming electron distributions in the left and right leads. Included in 
this self-consistent method is the constraint of particle conservation. We do not however 
include cumulative effects of the current on the phonon distribution in the conductor; i.e., 
we assume that the conductor is in contact with a heat bath. With the non-equilibrium 
electron distribution determined we then evaluate the current flowing through this system 
for finite temperatures and bias voltages. 

As an application of the above theoretical approach we calculate the electronic current as 
a function of voltage for a two site molecule. The molecule has two single-particle molecular 
energy levels through which electrons from the leads can resonantly tunnel. It also has two 
longitudinal phonon modes which can be excited. A SSH Hamiltonian is used to represent 
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the electron-phonon coupling. For our model system one of the phonon modes is chosen to 
have an energy equal to the gap between the two electronic energy levels. Because of this, 
electrons which are resonant with the higher energy level can emit a phonon and be resonant 
with the lower level. Such a process is found to be highly transmitting. When the voltage is 
high enough that there is a sufficient number of electrons in the left lead that can undergo 
this process, we find that the current in the right lead due to elastically scattered electrons 
is reduced significantly due to competition with these inelastically scattered electrons. Thus 
we find that the mutual exclusion between the scattered electrons that make up the non- 
equilibrium distribution in the drain lead can have important consequences for the electronic 
current fiowing through the molecule. 

In Sec. II we describe the class of systems that we study and explain how the scattering 
states are calculated. In Sec. Ill we present our generalization of the Landauer theory 
and our method for calculating the non-equilibrium the electron distribution and associated 
electronic currents. We then apply our methodology to a simple model for a two site molec- 
ular chain and calculate its current-voltage characteristics in Sec. IV. Our conclusions are 
summarized in Sec. V. 

II. INELASTIC SCATTERING IN ELECTRON TRANSPORT THROUGH A 

MESOSCOPIC CONDUCTOR 

In the absence of inelastic scattering and electron correlations, electron transport through 
a mesoscopic conductor can be described in terms of the probability T{E, V) that a single 
electron with energy E scatters through the conductor from the source to the drain at 
an applied bias voltage When phonons (or vibrational modes) of the conductor are 
considered, an electron entering from the source can suffer inelastic collisions by absorbing 
or creating phonons before entering the drain. Each of these processes can be described by 
its own transmission probability. If we assume that the source and drain are single mode 
leads which are free from phonons, we can characterize the above single electron scattering 
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process by a transmission probability T^^j^(E, E'). This describes an electron that enters 
from the left lead (L) with energy E, suffers inelastic collisions on the conductor which 
change the phonon population on the conductor from a to a', and then scatters into the 
right lead (R) with energy E'. Similarly, an electron incident from the left lead can suffer 
a collision and be reflected back into the left lead. This is characterized by a reflection 
probability -R^'" {E,E'). Elastic scattering processes are included with a' being equal to 
a. For all processes the total energy is conserved, so E' plus the net energy of all of the 
created and destroyed phonons is equal to E. We now describe a method for calculating 
these transmission and reflection probabilities for a mesoscopic structure. 

A schematic of our model for the mesoscopic conductor, which consists of a source 
(left) and a drain (right) lead attached to a one- dimensional conductor/wire is depicted 
in Fig. |1|. We assume that the leads are ideal and that they only have one electronic 
channel. Phonons are considered to be only on the conductor and not on the ideal leads. An 
electron flows through the wire under an applied bias V. We assume that the voltage drop 
only occurs at the interface between the wire and the leads with the left lead held at —V/2 
and the right lead held at +V/2. We solve the transport problem by solving Schrodinger's 
equation for the many-body wavefunction of the electron-phonon system from which we 
extract the transmission and reflection coefficients mentioned above. These will be used in 
the calculation of the non-equilibrium electron distributions and the current which follows 
in the next section. 

The Hamiltonian H of this system can be split up into the Hamiltonians for the isolated 
systems consisting of the left lead Hi, the right lead Hji, and the conductor Hq. It also 
includes the coupling between the wire and the left and right leads, Vl and Vr respectively. 
Thus we write, 

H = Hl + Hr + Hc + Vl + Vr (1) 

We treat each of the above Hamiltonians in the tight binding approximation. We model 
the wire using a set of orthogonal atomic orbitals Each atomic state has an energy Cj 
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with a corresponding creation operator bf which creates an electron in that state. In our 
Hamiltonian for the wire an electron can hop from orbital i to all nearest neighbor orbitals j, 
and this is governed by the hopping parameter t^ j for the equilibrium configuration of atoms. 
The wire also has a discrete set of vibrational modes, each characterized by a frequency cuk- 
The phonons on the wire are described by the state \a) = \{nk}), where {n^} is the set of 
mode occupation numbers. A phonon in mode k is created by the operator a'^. To first 
order in phonon creation/destruction operators the electron-phonon interaction is given by 
ji j, and this governs the process of an electron hopping from site i to site j and exciting or 
de-exciting a phonon in mode k. The Hamiltonian of the A^-site wire that we consider is, 

N 

He = Y.^ihth + J2tl{btb, + h.c.) + EE7y(a^ + ak)(b+bj + h.c.) +Y.hu;^4s.^ 

i=l i^j k i^j k 

The leads are modelled as a chain of atoms each with a single atomic orbital |n), where n 
labels the site. For the left lead n ranges from — oo to — 1. On the right lead n ranges from 
1 to oo. On the left lead an electron is created on site n by the operator c+. The site energy 
is given by = e + eV/2. Electrons can hop to nearest neighbor sites and this is controlled 
by the hopping parameter Similarly on the right lead an electron can exist on site n in 
an orbital with energy er — e — eV/2. The nearest neighbor hopping is controlled by the 
parameter P^. Since we assume that there are no phonons on the leads the Hamiltonians 
for the left and right leads are given by 

1 1 

Hl^ei E ^n^n + PL (Cn Cn-1 + h.c) 

n=— oo n=— oo 

oo oo 

Hr^SrY^ c:^Cn + pRj2iCnCn-l + h.c) 

n=l n=l 

The wire is coupled to the leads via Vl and Vr. We assume that only the sites on the leads 
directly adjacent to the wire are coupled to the wire. Electrons can hop from these lead 
sites to any site j on the wire and this is governed by w^ij for the left lead and Wij for the 
right lead. We take the electron-phonon interaction between the leads and the wire to first 
order. This allows for electrons from the leads to hop onto the wire and excite a phonon in 
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mode k. This electron-phonon interaction is given by r^i^^ for the left lead and T'lj for the 



right lead. Thus the coupling potentials are, 



N N 



Vl 



E^-i,,(c^i&, +h-c.) + EEr-i,j(ai[ + ak)(clibj +h.c.) 




N N 



Vr 



E t^i,i(c|6, + h.c.) + E E rtj(4 + ak)(c|bj + h.c.) 




We now determine the scattering wavefunction from Schroedinger's equation, = 
where a labels the initial phonon state on the conducting wire. The wavefunction 
describes the many body system consisting of a single electron and the phonons. As shown by 
Bonca and Trugman, the scattering process can be represented graphically.0 The power of 
this representation is in that the many-body problem can be viewed as a multi-channel single- 
electron scattering problem. The a' channel in either the left or the right leads corresponds 
to an electron propagating in the single electron mode of that lead with the wire being in 
its a' phonon state. A scattering process which changes the phonon state from \a) = |{^fc}) 
to \a') = \{n'j^}) can then be represented graphically as an electron incident from the a 
channel and then scattering into the a' channel. Fig. ^ shows a graphical representation of 
the scattering channels for a simple single site wire described by the above Hamiltonian. An 
electron incident in the a channel in the left lead can scatter elastically or inelastically into 
any of the outgoing channels in the left and right lead. Channels are connected horizontally 
by the electron hopping part of the Hamiltonian whereas they are connected vertically by 
the electron-phonon interaction. 

With this picture in mind, we now write explicit forms for the wavefunction on the left 
lead (L), right lead (R) and the wire conductor (C). We assume that the wire is initially in 
the phonon state |a) and that there is a rightward propagating electron in a Bloch mode 
Ene*"^n^) with energy E in the left lead. The electron can then be reflected elastically or 
inelastically into any of the other channels in the left lead and propagate in a new Bloch 
mode with energy E' in the left lead with the wire in a new phonon state \a'). Thus, the 
wavefunction on the left lead will consist of the initial state plus a sum over all possible 
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reflected states with some unknown coefficients ra',a- On the wire the wavefunction will be 
a linear combination of the product of the atomic orbitals with all of the different phonon 
states. On the right lead the electron can be transmitted into any of the outgoing channels 
and propagate in a Bloch state with an energy E', leaving the wire in a phonon state \a'). 
Thus the wavefunction on the right lead will be a sum over all transmitted states with 
coefficients t„' „ to be determined. Thus we have 



with 



e'^'y In) 



n=—oo 

N 



l^")c = EE<"U) 

a' 3=1 



(2) 
(3) 
(4) 



n=l 



1*^ 



R 



(5) 



The total energy of the system is conserved, so E = E' + Y.kiP'k ~ '<^k)hjJk-i where E' is the 
energy of the scattered electron, y"' corresponds to the reduced wavevector of an electron 
propagating with energy E' and is determined from the condition E' — OLLjR^'^^iLiR cos(|/"') 
depending on whether it is in the left or right lead. 

Inserting the above wavefunction into Schrodinger's equation E.\^°') — yields a 

system of linear equations which can be solved numerically for the r„/ q,, and the 
at each energy E. For our particular model based on single channel leads we then determine 
the transmission and reffection coefficients for each inelastic channel using, 



Tli^{E,E') 



-"R 



vl 



'^a' ,a I 



(6) 
(7) 



where vf^ and v'^ are the velocity of the electron in the a' channel of the left or right lead 
respectively. 
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We then solve the corresponding problem for electrons incident from the right lead to 
determine T^'_^j^(E, E') and -R^'" E'). We now proceed to show how these quantities are 
used to calculate the non-equilibrium electron distribution and current for the above system. 

III. NON-EQUILIBRIUM ELECTRON DISTRIBUTIONS AND THE 
EVALUATION OF THE ELECTRONIC CURRENT 

As mentioned in the introduction, when calculating the current using Landauer theory 
the Pauli exclusion principle for final states does not play a role; scattered electrons do not 
compete with each other. A complication that arises when inelastic scattering processes 
are considered is that electrons can now compete with other electrons that are scattered to 
the same final state; the Pauli principle becomes important. For instance, for the system 
discussed in Sec. II, an electron which suffers an inelastic collision that puts it into a final 
state with energy E' in the drain competes with all other scattered electrons that can occupy 
that final state (as shown in Fig. The likelihood that this state is unoccupied in the drain 
is (1 — /(£")), where f{E') is the non-equilibrium electron distribution function in the drain. 
Using the equilibrium electron distribution (Fermi function) here is not correct since we are 
discussing transport, an inherently non-equilibrium phenomenon. We will now explain how 
to determine this non-equilibrium distribution for the model presented in Sec. II. 

Fig. 1^ depicts the processes we are considering for electrons that scatter from various 
initial states with differing energies E to the same final state at energy E' in the right 
lead. Some assumptions must be made about which processes can compete with each other. 
To be consistent with Landauer theory we assume that elastic scattering processes do not 
compete with other elastic scattering processes (these correspond to those processes which 
both begin and end at energy E' labelled Al and Bl in Fig. ^). With this assumption an 
electron which is elastically scattered from the source to the drain at energy E' does not 
compete with an electron which is elastically reflected back into the drain at the same energy 
(process Bl). However, elastically scattered electrons can compete with all electrons that 
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can be inelastically scattered to the same final state. The inelastically scattered electrons 
correspond to those electrons which start at various energies E labelled A2, A3 and B2 in 
Fig. ^. We also assume that inelastically scattered electrons compete with all other scattered 
electrons that start in a different initial electron state (i.e. electrons undergoing process A2 
compete with electrons suffering processes Al, A3, Bl and B2). 

At a given final scattered energy E', we consider the following contributions to the 
non-equilibrium rightward moving electron distribution in the right lead: elastically scat- 
tered electrons transmitted from the left lead {f"''^i)^{E'), inelastically scattered electrons 
transmitted from the left lead (the final phonon state is a'), elastically scat- 

tered electrons which originate in the right lead and are reflected back into the right lead 
{f-ei)^{E')> inelastically scattered electrons which start in the right lead and are re- 
flected back into the right lead {f-^n)^{E'). (The (+) denotes transmitted electrons, the 
(— ) reflected electrons and the R signifies that these are the distributions in the right lead). 

To determine the contribution to the non-equilibrium distribution in the right lead that is 
due to transmitted electrons we must first consider the distribution of electrons incident from 
the left lead. As in Landauer theory, we assume that this distribution is given by the equi- 
librium electron distribution function which exists deep in the left lead. This is the Fermi 
function for the left lead El{E) = l/{exp{{E — ^L)/kT) + 1) where the electro-chemical 
potential is /iL = + eV/2 with ep being the common Fermi energy of the two leads and 
V the applied bias voltage. An electron incident with energy E from the left lead can then 
scatter elastically or inelastically into the final state with energy E' with a transmission 
probability T^!^j^{E, E'). (Elastic scattering corresponds to a = a' and E = E'). Because 
we now have a number of processes which can place electrons into this same final state the 
Pauli exclusion principle must be considered. To incorporate this into our determination of 
the distribution function we must then consider the probability that this state is unoccu- 
pied by competing electrons in the drain. For the distribution due to elastically scattered 
electrons, this probability will be determined by the non-equilibrium electron distributions 
of the inelastically transmitted and reflected electrons, according to the above assumptions. 

12 



It is given by (1 - Ea'(/+;»°n)^(^') - Ea'(/-;»°n)^(^'))- For the distribution of a specific in- 
elastic transmitting process, the probabihty that the state in the right lead is unoccupied by 
competing electrons is determined by the distributions due to all other scattering processes. 
It is given by (1 - {f^%nE') - {f^%nE') - EAf-^nE') - EpM^inE')). The 
transmitted electron distributions will be proportional to the product of the above probabil- 
ities with Fl and with the respective transmission probability T. We define the appropriate 
proportionahty constant for electrons incident from the left lead with energy E, to be c{E). 
This proportionality constant will be determined by imposing the constraint of particle con- 
servation. Combining the above considerations we arrive at the following expressions for the 
contributions to the non-equilibrium electron distribution due to elastically and inelastically 
transmitted electrons: 

imnE') = c{E')F,{E')n'^^{E' , E'){i - Y^uiinE') - e(/-£)^(£^o) (§) 

a' a' 

{&{E') = c{E)F^{E)T2^^^{E,E'){l - {f^nE') - {rjTiE') - 

Uf-'inE') - y: ininE')) (9) 

a' pj^a' 

Similarly for the contributions to the non-equilibrium electron distribution in the right lead 
due to elastically and inelastically reflected electrons, we must consider the distribution of 
incident electrons given by the Fermi function Fr{E) — l/{exp{{E — jiRj/kT) -\- 1) for 
the right lead where iir — ep — eV/2. Those electrons which are incident with energy 
E can be inelastically or elastically reflected into the flnal state with energy E' with a 
probability E!^" {E,E'). Again we take into account the competition for the occupation 
of this state using the Pauh principle. For elastically reflected electrons the probabihty 
that this state is unoccupied is the same as that for elastically transmitted electrons. For 
inelastically reflected electrons the probability that the state is unoccupied is determined 
by all distributions due to different scattering processes. It is given by (1 — (/"'g;)^(£") — 
(r;:j«(E')-E.'(/T£)''(^')-E/5^.'(/-;fn)''(i^'))- Again, the reflected electron distribution 
will be proportional to the product of the above probabilities with R and Fr. For electrons 
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incident from the right lead with energy E, the proportionahty constant is d{E). Based 
on these considerations we write the non-equihbrium distributions in the right lead due to 
elastically and inelastically scattered electrons as 

{tjT^E') = d{E')Fn{E')Rr{E\E'){l - Y^iftfufiE') - T^if-fufiE')) (10) 

a' a' 

if-inE') = d{E)Fn{E)Rr'{E,E'){l - UtT{E') - {t^^^E') - 

nuiinE') - if-inE')) (11) 

a' f3j^a' 

Similarly a system of equations is set up for the non-equilibrium electron distributions 
in the left lead, if"%)^, if+,in)^^ if-,ei)^^ and U-^)^- (Here, (+) signifies transmitted 
electrons from the right lead to left lead, (— ) refers to reflected electrons in the left lead, 
and (L) denotes that all these distributions are in the left lead). 

The last consideration we must impose is particle conservation. Particle conservation re- 
quires that the total incoming current be equal to the sum of the total transmitted current 
plus the total reflected current. This will allow us to determine the constants of proportion- 
ality for both the left and right leads. For the left lead we apply the following equality, 

fl{e) = ifz^E) + {f-%nE) + Y^utinE') + Y.u-if{E') (12) 

a' a' 

This equality expresses that the incident electron distribution is equal to the sum over all 
transmitted and reflected electron distributions that arise from that incident channel. Since 
all of the transmitted and reflected distributions originated from the same incident channel, 
they all involve the same proportionality constant c{E). Inserting the expressions for the 
/'s into Equation (|l^) allows us to solve for c{E). This yields, 

^ ^ Ael + Ain + Bel + Bin ^ ^ 

where 

Aei = n'ZAE, E){i - Yi&iE) - Yif-ifm 
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a' 

Bel = Rr{E,E){l - Y.UZn)'{E) - Y.U-i)\E)) 

a' a' 

Bin ^Y.Rr\E,E'){l - ifZnE') - ifZifiE') - 

a' 

Y.ir4n)\E') - E {r-i)\E')) 

Similarly, d[E) can be determined by applying a similar equality for the distributions origi- 
nating from the right lead. This yields, 

^ ^ Cel + Cin + Del + -Djn ^ ^ 



where 



Cei = T^'ZdE, E){i - Y.{r4^\E) - Y.u-t)'m 

a' a' 

Cin = Y.T^^l{E,E'){1 - {fl'%f{E') - ifZifiE') - 

a' 

Eif-D'iE') - E if+fyiE')) 
Del = RTiE, E){i - Y^UZufiE) - Y.U-ifm 

a' a' 

Din ^J:RT'{E,E'){1 - if^XiE') - it-TiE') - 

a' 

mninE') - e if-inE')) 

The above self-consistent system of non-linear equations must be solved iteratively to 
determine the unknown distributions /, and the unknown proportionality constants c{E) 
and d{E). For an initial trial solution, we let the /'s take on the following values, 

^ja,ajii^^-j _ El{E)T1'ZAE , E) 



TtCR{E,E)+Rr{E,E) 

(fOc,a-.R,T^-. ^ Fr{E)Rr°'{E,E) 

U-,ei) { ) T^'»^{E,E)+R'^"{E,E) 

/ .a,aNL/pN ^ Er{E)T^:^j;^{E , E) 

U+,ei) ^ ) T^:^^{E,E)+ir^''{E,E) 

(f»-»\L(p\ ^ Fl{E)RI''^{E,E) 
U-,ei) y ) T^;^j^{E,E)+Rl'''{E,E) 
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ifiinE) = 

{f^i)\E) = 

Using this initial guess for the /'s, we evaluate the constants of proportionality c{E) and 



d(E) from Equations |l^ and 14. With these determined, we then evaluate a new set of /'s 
from Equations |-|TI using these c's and d's and the old /'s. We then iterate until the /'s 
have converged. 

With the non-equilibrium electron distributions thus determined we proceed to evaluate 
the electronic current flowing through the wire for an applied voltage V and at a given 
temperature T. At a given energy E the current, {6i+)^, flowing in an energy interval 5E 
due to rightward propagating transmitted electrons in the right lead is given by the charge 
— e times the velocity vr{E), times the density of states dg{E)/dE, times the electron 
distribution for transmitted rightward propagating electrons in the right lead {f+^tot)^{E) = 
(/+'e«)^(-^) + J2a'if+An)^i^)- Thus wc cau Write (including a factor of 2 for spin) 

i6t^)'' = -^iU,totnE)6E (15) 

where in one dimension the density of states cancels the velocity factor. The total current 
due to rightward propagating transmitted electrons in the right lead is found by integrating 
Equation yielding, 

(h)^(V^) = -jJdE {U,totf{E) (16) 

Similarly, the total current in the left lead due to transmitted leftward propagating electrons 
is given by, 

innV) = jJdE {U,totf{E) (17) 

where [f+^totf{E) = {f+%)^{E) + Ea'{f+fn)^{E). The net current flowing through the 
mesoscopic conductor is then 

ttotiV) = i^^nV) + i^^nV) (18) 
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IV. APPLICATION TO A DIATOMIC MOLECULE 



We now proceed to apply the above formalism to a simple molecular wire system. It 
consists of two identical ID leads which are attached to the left and right ends of a two 
atom chain. Fig. ^ shows a schematic of our system. Each atom making up the molecule 
has just one atomic orbital, so the isolated two atom molecule has two discrete electronic 
energy levels. When this molecule is coupled to the leads, incident electrons with ener- 
gies that roughly coincide with the energy of either of these two states will be favourably 
transmitted. In the absence of phonons we expect two strong peaks in the transmission 
spectrum at the energies of the two molecular levels. When longitudinal phonons on the 
molecule are considered we expect interesting physics to occur when one of the phonon en- 
ergies corresponds to the energy gap between the two energy levels of the molecule. When 
an electron is incident at an energy corresponding to the higher molecular level it can emit a 
phonon and drop to the lower electronic energy level. This process should have a significant 
transmission probability because the electron will be resonant with the higher energy level 
and then also with the lower energy level. (Processes where phonons are created and the 
electron is not resonant with either of the two molecular energy levels will not be strongly 
transmitting). If we assume that the leads' Fermi energy lies between the two molecular en- 
ergy levels, we expect the following behavior for the non-equilibrium electron distributions: 
For low voltages and temperatures, the non-equilibrium electron distribution in the drain 
will be dominated by resonant electrons elastically scattered through the lower energy level 
and elastically reflected electrons. As the voltage increases the electro-chemical potential 
in the left lead rises while it decreases in the right lead. Eventually, at sufficiently high 
voltages there will be electrons in the left lead with enough energy to undergo the inelastic 
scattering process mentioned above. The non-equilibrium electron distribution in the drain 
will then have a significant contribution due to these inelastically scattered electrons. These 
inelastically scattered electrons will compete with elastically scattered electrons that are 
transmitted through the lower energy state. This competition between electrons due to the 
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Pauli principle will play an important role in determining the current of this two site system 
at higher voltages. 

We have chosen the following set of tight-binding parameters for the system described 
above. On the leads, the site energy is e = —10 eV and the nearest neighbor hopping 
energy is (3l = I3r = —2 eV. For the molecule, the two atoms have the same site energy 
ei = 62 = —10.0 eV and the hopping energy is ti^2 = —0.15 eV. The coupling of the 
left lead to the first molecular site is = —0.2 eV and the coupling of the second 

molecular site to the right lead is W2,i = —0.2 eV. We use the SSH Hamiltonian to model 
the vibrational modes of the molecule and electron-phonon coupling as is discussed in the 
Appendix: Our vibrational model of the molecule has two longitudinal phonon modes with 
energies hui = 0.24 eV and huj2 = 0.33 eV. The electron-phonon (e-p) couplings on the 
chain are 7i2 = 0.0 eV and 7i2 = —0.24 eV. The e-p couplings between the leads and the 
chain are = 0.12 eV, = 0.12 eV, ^ = -0.12 eV and ^ = 0.12 eV. We assume 
that the system is at a temperature of T = 77 K. 

For the transport calculation, the initial phonon state \a) = \{n^}) is taken to be the 
average phonon distribution on the chain for temperature T, with n'' = l/{exp{hujk/ {ksT)) — 
1). We approximate each average mode occupation number by the nearest integer. At 
the temperature of interest, \a) = \{ni,n2}) = |{0,0}). The other phonon states that we 
include in the calculation are those obtained by taking all rij = 0, . . . , N^ax where N^ax = 2. 

In Fig. ^a, the total transmission probability found by summing over all outgoing channels 
(solid line) is shown along with the transmission probability in the absence of phonons 
(dashed line). (These transmission probabilities were calculated for V applied across the 
molecule. We emphasize that they are the one- electron transmission probabilities defined 
in Section 11 and do not reflect the effects of the Pauli principle discussed in Section 111. 
The latter is taken into account in the calculated non-equilibrium distributions and currents 
presented later in this Section). In the absence of phonons the two strong resonant tunnelling 
peaks are seen. With the inclusion of phonons, there are extra peaks due to the creation of 
phonons. The two resonant peaks are shifted because the inclusion of vibrations alters the 
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hopping between sites which in turn shifts the energies at which electrons are resonant. In 
Fig. ^ the transmission probabihty for elastically scattered electrons is shown. The elastic 
spectrum show two strong peaks at the molecular energy levels. It also has other smaller 
peaks which are due to the interaction of the electrons with virtual phonons. 

Fig. 1^ shows the inelastic one-electron transmission probability for three of the inelastic 
channels vs. the energy E of the incident electron state. The first graph, Fig. ^a is for 
the process where a phonon with hui = 0.24 eV is created on the molecule. The peak 
labelled A corresponds to an electron emitting this phonon and dropping into the lower 
resonant molecular energy level at -10.26 eV. The peak labelled B corresponds to resonant 
tunnelling through the higher molecular energy level after the creation of this phonon. The 
two tunnelling peaks C and D at the energies of the molecular resonant states correspond 
to an electron entering the resonant molecular energy levels prior to emitting the same 
phonon. The second graph, Fig. ^ corresponds to the creation of a single excitation of 
the second phonon mode with energy huj2 = 0.33 eV, equal to the molecular energy level 
gap. As discussed previously, we expect electrons with energies incident on the higher energy 
molecular state to be strongly transmitting and this process corresponds to the peak labelled 
A in this graph. Graph (c) is for the process where a single excitation is created in both 
phonon modes. The transmission probability of this process is essentially the superposition 
of Fig. ^ and b, with an overall reduction in amplitude because this corresponds to a higher 
order process. 

We now proceed to calculate the non-equilibrium electron distributions from the cal- 
culated one-electron transmission and reflection coefficients using the method presented in 
Sec. III. We take the Fermi energy of the leads to be ej? = — 10 eV. The non-equilibrium 
distributions play an important role in the calculation of the current because of the Pauli 
principle. For instance, the presence of inelastically transmitted or back-scattered electrons 
can result in a significant suppression of the elastically transmitted current. In this sys- 
tem, as mentioned above, the process whereby an electron creates a phonon with energy 
huj2 and is then at resonance with the lower electronic resonant tunnelling state, is strongly 

19 



transmitting in the inelastic spectrum at E = — 9.9eV. At low voltages there are not very 
many of these particular inelastic processes occurring since there are not many electrons 
with this energy incident from the left lead, and the final electron states for the transitions 
are almost fully occupied by competing elastic processes. At higher voltages (greater than 
0.55 V) there is a significant number of these processes and the non-equilibrium distribution 
f1'^^{E) where a' corresponds to the phonon state {^^} = {0, 1} is significant. As the volt- 
age increases further, the total distribution due to inelastically scattered electrons continues 
to grow as more processes become energetically favourable. The total non-equilibrium Fermi 
distribution for all of the transmitted elastically scattered electrons is shown in Fig. |^a for 
0.5 V (solid fine) and 1 V (dashed line). Fig. 0b shows the corresponding distribution for 
inelastically transmitted electrons. It highlights how this distribution function increases as 
more inelastic processes become energetically favourable. It also shows that there is a sig- 
nificant presence of inelastically scattered electrons in the drain for energies around -10.25 
eV. At this energy, elastically scattered electrons are also particularly strongly transmitting 
and the dip in the distribution of the inelastically scattered electrons at this energy is due 
to the presence of these elastically scattered electrons in the drain. This is labelled "A" in 



These non-equilibrium distributions are now used to evaluate the current fiowing through 
the two atom molecule. In Fig. Ha,b we show the total rightward fiowing current that is 
elastically and inelastically transmitted through the molecule, evaluated using our method. 

We also show (dotted lines) the corresponding currents evaluated assuming instead that 
the non-equilibrium distribution of electrons in the drain that appears in the Pauli exclusion 
factor is replaced simply by the equilibrium Fermi function in the right lead. For transmitted 
elastically scattered electrons this simplified expression for the current is 



For models that do not include electron-phonon interactions, the net current through the 
molecule obtained in this way is the sam as is obtained from Landauer theory, although 



Fig. 0b. 




(19) 
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the factor (1 — Fii{E)) does not appear in the latter as was explained in the introduction. 
However if electro n-phonon interactions are present, approximating the effects of the Pauli 
principle in this way does not account for competition between elastic and inelastic processes 
which (at higher voltages) inject electrons into states that are unoccupied according to the 
equilibrium distribution. This approximation (even at low voltages) also violates particle 
conservation when electron-phonon interactions are present. 

By using the non-equilibrium distributions described in Section III these effects are prop- 
erly taken into account and current is conserved. All of this is reflected in the graphs: In 
Fig. Ha initially the solid curve (the present result) shows significantly more elastic current 
(as expected from Landauer theory), and then eventually the elastic current is reduced due 
to the presence of transmitted inelastically scattered electrons in the normally unoccupied 
states of the drain. The total current flowing through the molecule as a function of voltage 
is shown in Fig. |c. The larger total current exhibited by the solid curve is a manifestation 
of current conservation in the present formalism. The dashed curve falls below the solid 
curve because the approximations used to calculate it do not conserve current. 

V. CONCLUSIONS 

We have shown that to generalize Landauer theory to include the effects of inelastic 
scattering, it is necessary to calculate self-consistently the non-equilibrium contributions to 
the electron distribution in the source and the drain leads. We have developed a method 
for determining these non-equilibrium distributions and applied it to a simple model. The 
model was based on single- channel ideal leads attached to a mesoscopic conductor. Phonons 
were only considered on the conductor. The non-equilibrium distributions were determined 
solely in terms of the Fermi functions of the left and right leads and the transmission and 
reflection probabilities for the various one-electron elastic and inelastic scattering processes 
considered. We showed how to properly take account of the Pauli exclusion principle and 
particle conservation in the calculation of the non-equilibrium electron distribution. We 
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applied our approach to a two site molecule. One of the phonons had an energy that 
corresponded to the molecular energy level gap. For this system we found that the details 
of the non-equilibrium electron distribution had an important role in the determination 
of the current. At higher voltages we found that the distribution in the right lead due 
to inelastically scattered electrons began to compete significantly with elastically scattered 
electrons because of the Pauli exclusion principle. We found that this competition resulted 
in a reduction of the elastically scattered current at higher voltages. 
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VII. APPENDIX 

To find the vibrational modes of a molecule the atoms and bonds may be modelled by a set 
of balls and springs. For our model, we consider chain-like molecules, and represent the chain 
as a set of balls with effective springs connecting nearest neighbor atoms. The longitudinal 
displacement from equilibrium of the i*'* atom with mass rrii in the molecule is given by 
the coordinate Qi. We assume that each atom only interacts with its nearest neighbours 
via an effective spring coupling Vij. The molecule has a discrete set of longitudinal normal 
modes which vibrate with a particular frequency Uk and have atomic displacements given 
by = exp {—icukt)- These normal modes satisfy the familiar eigenvalue problem, 

J2V^,,d'J = ulmd^ (20) 
j 

where for our model, the sum over j is just over the nearest neighbors of atom i. Solving this 
equation yields the eigen-frequencies and eigen- vectors for the longitudinal normal modes. 

The problem is canonically quantized by assuming [qi,Pj] = ihSij, where pj is the con- 
jugate momentum to qj. We transform into normal mode operators Qk and Pk by writing 
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Qi = J2kdiQk and pi = J^kdiTriiPk. It is easily verified that [Qk,Pk'] = ifiSk,k'- Performing 
second quantization on tliese operators gives, 



Qk = J^{ak + at) (21) 



Pk = ^\|^{at-a,) (22) 

where and are creation and destruction operators for the normal mode (or phonon) 
with frequency Uk respectively. 

The effects of longitudinal vibrations on the electronic properties of the molecule can 
be taken into account by using the Su-Schreiffer-Hieger (SSH) mo delJli In it, the nearest 
neighbor hopping parameter tij between two neighboring atomic sites is expanded to linear 
order in terms of their atomic displacements from equilibrium. This gives, 

where t^^ - is the hopping parameter at q = qi — qj = and = {dtij/dq)q=Q. (Note 

By using the above expression for qi in terms of the QkS, and writing them in their 
second quantized form, this gives tij as. 



kj = th - - 4)y + 4) (24) 



This defines the mode dependent electron-phonon coupling, as 



in terms of the known quantities aij, uj^ and the (i'^'s. 

We now proceed to evaluate the eigen-frequencies and eigen-vectors of the longitudinal 
modes for our two atom molecule attached to two semi-infinite leads. The atoms in the 
leads are assumed fixed in their static positions. The effective spring coupling between the 
leads and the nearest neighbor atoms on the molecule is given by = ISQeV/A^. The 
effective spring coupling between the two atoms of the molecule is = 90eV/A^. Each 
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atom of the molecule has a mass m = 13 A.M.U.. Solving the eigenvalue problem using the 
above parameters we find two phonon modes with energies of hcui — 0.24 eV and hu!2 — 0.33 
eV. The corresponding eigen-vectors are d\ = 0.054^4, dl = 0.054A and = 0.054A, 
dl = -0.054A. 

In the SSH model the electrons are coupled to the phonons via the parameter aij. The 
nearest neighbor coupling between the leads and the molecule and between the two sites on 
the molecule is assumed to be — 25 eV/Afor i to the right of j. This gives the couplings 
between the first site on the molecule and the left lead as = T? -i = 0.12. The couplings 
between the right lead and the second site on the molecule are 2 = —0.12 and 2 = 0.12. 
The electron-phonon couphngs on the molecule are 72 1 = and ^21 — —0.24. 
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FIGURES 

FIG. 1. A schematic of our model for a one dimensional mesoscopic conductor. The conductor 
bridges two one dimensional ideal leads. The applied bias voltage is V. The Hamiltonians of the 
individual elements making up the system are labelled. 

FIG. 2. A graphical representation of the inelastic scattering problem for a single site wire with 
first order electron-phonon coupling between the leads and the wire . Each phonon state of the 
wire along with the Bloch state of the electron in a lead can be visualized as a separate scattering 
channel. Scattering channels are connected horizontally by the electron hopping parameters t and 
P, while they are connected vertically by the electron-phonon interaction F, 7. An electron is 
initially incident in the a channel and can scatter into any of the other channels. 

FIG. 3. A schematic of the various elastic and inelastic scattering processes that can occur 
within the mesoscopic conductor that contribute to the rightward propagating electron distribution 
in the right lead. The applied bias voltage is V. The electrochemical potentials of the two leads 
under this applied bias are labelled fiL and fiR for the left and right leads respectively. They share 
a common Fermi energy labelled by e^. Scattering processes originating in the left lead are labelled 
A1-A3, and for those originating in the right lead, they are labelled B1-B2. 

FIG. 4. A schematic of our two-site molecular wire. We model the longitudinal vibrations of 
the molecule using a nearest neighbor ball and spring model. The two atoms of the molecule are 

coupled by an effective spring coupling k, while they are coupled to the leads with a spring coupling 
K. 

FIG. 5. (a) Total transmission probability as a function of energy E of the incident electron for 
the two site molecule. The dashed curve corresponds to the transmission probability in the absence 
of phonons. The solid curve is with the inclusion of longitudinal phonons in the calculation at a 
temperature of T = 77 K. (b) The transmission probability for the elastic scattering channel of the 
two site. Both graphs are dX V = 0. molecule. 
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FIG. 6. Transmission probabilities for three of the inelastic scattering channels of the two site 
molecule, (a) Scattering from the initial state |q) = |0,0) to |a') = |1,0). (b) Scattering from 
the initial state |a) = |0,0) to \a') = |0, 1). (c) Scattering from the initial state |a) = |0, 0) to 
\a') = |1, 1). All graphs were calculated at T = 77 K and V = Q. 

FIG. 7. Non-equilibrium electron distributions in the right lead for the two site molecule, 
(a) Total non-equilibrium electron distribution for transmitted elastically scattered electrons at 
V (dashed curve) and 1 V (solid curve), (b) Total non-equilibrium electron distribution for 
transmitted inelastically scattered electrons at 0.5 V (solid curve) and 1 V (dashed curve). The 
postion labelled "A" shows how there is a dip in the inelastic distribution due to the sharp rise in 
the elastic distribution at that energy E. 

FIG. 8. (a) Total rightward flowing electronic current in the right lead due to elastically trans- 
mitted electrons, (b) Total rightward flowing electronic current in the right lead due to inelastically 
transmitted electrons, (c) The net electronic current flowing through the two site molecule. The 
solid curves are calculated from the non-equilibrium electron distribution and the dashed curves 



are calculated using Equation 1£. The current is in units of 2e /heV . 
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